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We analyze a two dimensional model of gauged fermions with quartic couplings in the large-TV 
limit. This combines the 't Hooft model and the Gross-Neveu model where the coupling constants of 
both theories are arbitrary. Analytic equations describing the meson states of the theory are derived 
and are solved systematically using various methods. The physics of the model is investigated. 
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I. INTRODUCTION 



Solvable models have greatly contributed to our understanding of the dynamics of quantum field theories. Two such 
1 + 1)- dimensional models solvable in the large N limit which are "classics" in this regard, are the 't Hooft model 
gauge theory with fermionic matter, and the Gross-Neveu model [^|, a model with a four fermi coupling. (For 
reviews on the subject, see, for instance, S.) Both these models are "solvable" from first principles yet they are not 
I . integrable models in the usual sense, except for the case of massless Gross-Neveu model [ 3 . These types of models 
are hard to come by and we believe that they contribute to the understanding of more realistic theories. Indeed, the 
't Hooft model and the Gross-Neveu model have a wide area of applicability, as evidenced by the contribution of these 
models in many areas of physics, including but not restricted to particle and nuclear physics, even recently. In this 
paper, we solve a model in the large-iV limit wherein the four fermi couplings and the gauge coupling coexist with 
£f~) ' arbitrary strengths, thereby enlarging this class of models in an essential way. We extend and generalize the work of 
Burkardt ||] , who derived the meson bound state equation in the gauged four fermi model and analyzed the equation 
. from a somewhat different approach from ours. 

The contents of this paper are as follows: First, in section 2, we solve the Gross-Neveu model using the Bethe- 
Salpeter equation, which, to our knowledge, has not been presented previously. In this section, we fix the notation and 
summarize the physics of the Gross-Neveu model, which will be useful later on. In section 3, we analyze the theory 
Oh! of gauged fermions with four fermi couplings in the large-TV limit. The model combines and generalizes the model 
of 't Hooft and the model of Gross and Neveu. The model is more general than the Gross-Neveu model even when 
\ the gauge coupling is zero. Next we derive analytic equations for the meson states of the theory. Renormalization 
i prescription is given and equations involving only finite physical parameters are derived. In section 4, methods are 
presented in detail for solving the meson state equations systematically. Using these methods, we derive various 
results on the physical properties of the model which are analyzed in section 5. We end with a summary and a more 
general discussion regarding the model in Section 6. 



II. THE MASSIVE GROSS-NEVEU MODEL 

In this section, we analyze the Gross-Neveu model and derive the Bethe-Salpeter equation for the fermion- 
antifermion (meson) channel in a spirit similar to that of 't Hooft's analysis of two dimensional QCD fjl]. This 
method is distinct from the methods applied to the Gross-Neveu model previously and will provide useful par- 

allels for some aspects of the model to be discussed below. We compare the results obtained from the Bethe-Salpeter 
equations to those obtained from the usual approach and also summarize the physics of the model pertinent for the 
sequel. 

The Gross-Neveu model we analyze in this section has the following Lagrangian; 
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This is equivalent to the following Lagrangian written using the auxiliary real scalar fields a, x, 
L = jfifaj - - {a 2 + x 2 ) + aW {a + ixj 5 ) ~ ~<r j = 1,2 . . . , N 



(2) 



We will use light cone coordinates v 
become simple in the chiral basis 



A. The equation for the "meson" bound states 

v T = 1/V2 (y° ± v 1 ) below. In light cone coordinates, the gamma matrices 
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l/N expansion is performed by expanding in powers of 1/N while treating a 2 N as a quantity of order one. We use 
the metric (+, — ). 

From the interactions in the Lagrangian (^) we obtain a self consistent equation for the propagator in the large 
N limit. As is clear from the interactions, these corrections do not have any momentum dependence. Therefore, it 
can be effectively summarized in a mass parameter, which we shall call M. The fermions in this model are physical 
particles and the parameter M is the physical mass of these fermions, which will be determined self consistently. 
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FIG. 1. Bound state equation for the fermion-antifermion system in the Gross-Neveu model 



Using the full propagator, we may straightforwardly obtain the Bethe-Salpeter equation for what we shall call the 
"meson wavefunction" , ip a p, from the contributions graphically represented in fig. [J as: 



^af){p,r) = 2a 2 b N 
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S(p) is the full fermion propagator in this model, which is none other than the tree-level propagator with the mass m 
replaced by M. a& denotes the bare four fermi coupling. Following 't Hooft, we integrate over one component of the 
momentum p + and define (p a f3 = J dp+ipap. Then, after some computation, the equations for the meson wavefunction 
ij) a p reduce to the following equations for the components: 



2tt 



2tt 



2tt M 
ajNy/2r- 

2tt V?r- 
ajN M 



P 2 



A 2 



A 2 



P 2 



1 



x(l — x) 



(pix{x) = - 



1 



x(l — x) 



1 

x(l — x) 

1 

x(l — x) 



2(1 -x) 



A 2 



l-x 



^22(2;) = - — 

2x 



012 (x) 



x ( l - x) J 

/1 2 + - 

X 



1 



(i 



x(l — x) 



l-x 
1 

dy<pu(y) 



dy (P22 (y) 







1 - X 



dy<pn{y) 



1 



x 



dy$22{y) 







1 



2x(l - x) 



~2 1 1 

X 1 — X 



dy<f>u{y) 



2x(l - x) 



1 

X 



l-x 



dy ^22 (y) 



(5) 



2 



Here, we defined the momentum fraction of the incoming antifermion x = p~/r_ and the mass squared of the bound 
state in units of the fermion mass squared as /2 2 = 2r + r_/M 2 . Without any loss of generality, we may define 
j dxtfn = 1, j dx(f>22 — C, where C is to be determined later. Then all the meson wavefunctions (p a p may be solved 
algebraically using the equations (^|) as follows. The consistency with the normalization of (pn requires that 

4tt f 1 - - 2 + o-7T \ 

dx—, .„ 2x(1 ~ x) (6) 



a 2 N J 1 - fl 2 x{l - x) 

The compatibility of this with the normalization condition for <f2i requires that C = ^fl. These two cases correspond 
to the meson states \ an d cr respectively. We obtain the equations determining the masses of x an d a as 

X ' a 2 N J o aX l-p? x x{\-xy °- a 2 N J, 1 - ftx(l - x) ( '> 
Here, we renormalized the coupling constant as 

4tt 4tt 1 f 1 1 

The integral needs to be regularized at the endpoints 0, 1 which is not explicitly expressed here. The same renormal- 
ization was employed in the light front Hamiltonian approach in [gj . This regularization is effectively an ultraviolet 
cutoff, which will become clear below. The meson wavefunctions for o can also be obtained algebraically as 

1 1 — 2x 

X- ^ 2 (i)=const.X- ^— -, a: ^ 2 (x) = const, x — (9) 

1 — fi^xyl — x) 1 — fL„x\l — x) 

The component if 12 is shown here since it corresponds to the relevant component of the meson wavefunction in 
the 't Hooft model |Q and will also be the essential component in our analysis of the gauged four fermi model. 
The wavefunction for \ IS consistent with the results obtained using light front Hamiltonian methods Other 
components of the wavefunction can also be computed algebraically using (^|). 



B. The analysis of the Gross— Neveu model using scalars and its relation to the Bethe— Salpeter equation 



In this section, we clarify the relation between the results obtained above using the Bethe-Salpeter equation and 
the results obtained from using the auxiliary scalar fields cr, x as m t ne Lagrangian (^|). Since the approach of using 
scalars is more standard and is explained elsewhere, we refer the derivation to the literature (|J^|. 

The full propagators for the a, x fields are 
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D- 1 (p 2 ) = l 
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where the a& is the bare coupling and A is the ultraviolet momentum scale cutoff. (There is a mild abuse of notation 
here; this bare coupling is in principle not the same as the one in the previous section since the regularization methods 
are different.) The function B(p 2 ,M 2 ) is defined as 



B(p 2 ,M 2 ) = y/l -4M 2 /p 2 In 
The renormalized coupling a defined in (& 



V / l-4M 2 /p 2 + l 
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y/1 -AM 2 /p 2 - 1 
is related to the bare coupling ai, here as 
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We find that the equations for the poles in the propagators for cr, x in ([h]) indeed agree with the equations (0) in 
this renormalization scheme. The propagators have cuts for p 2 > AM 2 corresponding to the production of physical 
fermion— antifermion pair of mass M each. 
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The effective potential for the scalars may be obtained by computing the contributions from the fermion loops in 
the large- TV limit as 

1/9 o\ m a?N , , os ( ai (a 2 + y 2 ) \ „ 
V{a, X) = - 2 (^ + X 2 ) - -a + -A- (. 2 + X 2 ) (in - ij (13) 

By minimizing the potential, we obtain a vacuum expectation value (a) for cr. The physical mass is related to the 
vacuum expectation value simply as M = a&(er), without loss of generality. We obtain a relation between the bare 
and the renormalized parameters 

Mm 

Since the mass is dynamically generated in the Gross-Neveu model even when m = 0, the chiral limit corresponds to 
a 2 N -> oo. 



C. Physics of the Gross-Neveu model 



Here, we briefly summarize the physics of the Gross-Neveu model, in particular, emphasizing the salient features 
and its underlying physics which will be useful later on. The Lagrangian (Q) has two parameters, m and a. Due to 
dimensional transmutation, the model is determined essentially by only one parameter. We can solve the equations 
(0) or the pole equations of the propagators (|o|) to obtain the masses of the scalars a, \- We plot the spectrum of 
cr, x against the inverse coupling in fig. 




We understand the various regions in the coupling constant as follows: 

1. l/a 2 N — 0: The chiral point. Here, the masses for a and x are respectively 2M and zero. The wavefunction 
for x, (p*2 [ x )i is a constant in this limit, similarly to the 't Hooft model, x is the "Nambu-Goldstone" boson 
of the theory. Strictly speaking, Nambu-Goldstone boson does not exist in (1 + 1) dimensions ||, yet it is 
well known that many physical aspects of the higher dimensional theories are also seen in (1 + 1) dimensional 
theories, especially in the large N limit. A similar massless boson bound state exists in the 't Hooft model in 
the chiral limit. 
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The status of the a particle is interesting; while the a particle is usually said to exist, its wavefunction <fi 2 (x) 
approaches const. /(I — 2x) and is singular in the limit l/a 2 N — > 0. Physical decay into a fermion and an 
antifermion becomes kinematically allowed in this limit /j 2 — » AM 2 and the singular behavior is due to this. 
Similar behavior is also seen for \ m the limit fj, 2 — > 4M 2 . The wavefunction is well behaved when a 2 N < 0, 
yet in this region, the vacuum is not physically stable, as explained below. 

2. a 2 N > 0: This region is physically consistent. The mass of x is between zero and 2M. The meson wavefunction 
has a singular limit (p\ 2 — >const.x(l — 1x)~ 2 as l/a 2 N — > oo. The pole in the a propagator (|t]) that exists for 
a 2 N < ceases to exist in this regime and there is no bound state corresponding to u. Poles corresponding to 
resonance states also do not exist, unlike the Gross-Neveu model in (2 + l)-dimensions ||. 

3. a 2 N < —ir: While a scalar has a finite mass, x is tachyonic. From the potential, we may understand this as 
follows; we are at an unstable vacuum where the potential is locally a minimum in the a direction yet maximal in 
the x direction. Choosing the physically sensible vacuum reduces this case to the previous physically consistent 
case. 

4. — 7r < a 2 N < 0: We again have chosen an unstable vacuum, this vacuum is unstable in both a and x directions 
so that both a, x are tachyonic. Had the correct vacuum been chosen, the theory reduces to the a 2 N > case 
above. 

We plot the potential for these various cases in fig. [1 along the plane x = 0- 




example is given for the each of the four cases 1 — 4 explained in the text. The values chosen are (2n)/a 2 N = 0,1,-1 and -3, 
respectively. Vertical scale is in units of NM 2 /2%. 
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III. THE GAUGED FOUR FERMI MODEL 



In this section, we analyze the gauged four fermi model. We obtain the Bethe-Salpeter equation for the fermion- 
antifermion bound states — which we shall call the "mesons" for obvious reasons — and perform the renormalization 
of the model. We further systematically solve the Bethe-Salpeter equation to obtain the masses and the wavefunctions 
of the meson states. 

The Lagrangian of the gauged four fermi model we work with is 

C= -~tr (F^Fn + Yj>f{0-™f)i>f + y E {W>^fWfi>f>) " y E ^hl^fWn^f) (15) 
/ (/,/') (/,/') 

The covariant derivative is defined as = — igA^ 1 where g is the gauge coupling constant. The color indices 
have been suppressed. Index / denotes the "flavor" index and is included here since we shall consider bound states 
involving fermions of different masses. The motivations for such a generalization is obvious when we want to apply the 
model to more realistic situations. The Lagrangian generically respects the global flavor symmetry [U(l)v] N , which 
is enlarged to SU(N)v X U(l)v when all the fermion masses to/ are the same. This symmetry is further enlarged to 
the chiral flavor symmetry SU(N)l x U(1)l x SU(N)r x U(1)r when all to/ = 0. The Gross-Neveu model we dealt 
with in the previous section corresponds to the case when there is no gauge coupling, only a single flavor type and 

2 2 



A. The equations describing the mesons 

Below, we will derive Bcthc-Salpctcr equations for the fermion bound states by using methods similar to those of 
't Hooft Q. The situation is intrinsically more complicated than that of 't Hooft since the four fermi interactions 
mix the fermions of different chirality and the Bcthc-Salpctcr equation can no longer be straightforwardly reduced 
to a one dimensional equation. We fix the gauge to the light-cone gauge A_ = A + = 0. Light-cone gauge has the 
advantage that there are no gluon self-interactions in (1 + 1) dimensions. We take the large N limit by letting N go 
to infinity keeping g 2 N, a 2 N, a 2 N fixed. 

First, we obtain the full propagator self consistently from the Schwinger-Dyson equation. Diagrammatically, the 
equation may be expressed as in fig. in the large TV limit. 
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FIG. 4. The self consistent equation for the propagator in the gauged four fermi model 



Solving the equation we obtain the full propagator S(p; Mf) as 

g 2 N ( sgn{pJ) 



S-\p\M f ) 
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where Mf is the mass parameter containing the quantum corrections, as in the previous section. We introduced a 
cutoff A_ for the infrared divergence in p_ integral. 

The bound state equation for fermion- antifcrmion bound state may be derived in the large iV limit, which is 
diagrammatically depicted in fig. |^. 
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FIG. 5. The bound state equation for the fermion-antifermion system in the gauged four fermi model. 
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Denoting the wavefunction of the bound state as ^ a/ 3(p,r), the Bethe-Salpeter equation can be derived as the 
following matrix equation in the large N limit 

^(p,r) = -i2g 2 NS(p;M 1 ) 1 + JAL ^(p + k,r)j+ S(p - r; M 2 )-^ 
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tr (j 5 tp(k,r)) 



(17) 



The suffix b on the couplings indicates that these couplings are bare parameters. We define, as in the previous section, 
<Pap{P-,r)= dp+ip a p(p, r). Then the bound state equations may be obtained for the components as 
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/^-^ M 2 --3^' *=^> ( 22 ) 



Here we defined 

M 2 2 = 7r2r + r - = ^ 
V^' M ~ <? 2 iV ' X " r_ 

and Pj denotes the principal part integral. When the four fermi couplings a 2 , a 2 , are absent, the equation (J2TJ) is 
the 't Hooft equation, which is a closed equation by itself. Here, all the equations are coupled and they need to be 
disentangled in a more sophisticated manner. 

Superficially, we have as many equation as the unknowns — the meson wavef unctions, <p a /3's. However, we expect 
the infrared cutoff A_ to decouple from these physical equations, so that these equations are possibly over-constrained. 
It may be shown that all these equations consistently reduce to the following equations ( |23| ) — (H|) and (^0|). These 
equations do not involve the infrared cutoff but are yet to be renormalized: 
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From these equations we may derive the following closed equation for ip\2 (= tp), whose suffix we shall omit for brevity 
from now on. 

H 2 <p(x) = H(p(x) 
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We shall often refer to H as the Hamiltonian. It is clear that this equation reduces to the 't Hooft equation when 
a 2 = a 2 = and that it reduces to the Gross-Neveu model case obtained in the previous section when g = 0, fa = fa 
and a 2 = a 2 . Even when the gauge coupling is zero, this model is more general than the massive Gross-Neveu model 
in that it incorporates different four fermi couplings a, 0,5 and fermions of different masses. This equation describes 
the properties of the mesons in the gauged four fermi model. When the fermion masses are equal, fa = fa(=fa), the 
equation takes a substantially simpler form; 

V = ~Tx \V\ X ) ~ P I d V 
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(27) 



When the fermion masses are equal, fa = fa, and when the couplings are equal, a 2 — a 2 , the bound state equation 
for the mesons in the gauged four fermi model (p6| ) reduces to the equation given by Burkardt || . Burkardt obtained 
a renormalized form of the equation when the meson wavefunction is an even function of the momentum fraction by 
using an operator identity involving the divergence of the axial current. In contrast, below, we will renormalize the 
more general bound state equation and reduce the equations to its renormalized form without using any further 
identities. 
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B. Renormalization of the gauged four fermi model 



The equations derived in the previous section (|23|)-(|26|) are expressed in terms of bare quantities. The equation 
for the meson wavefunction should be expressible in terms of renormalized quantities and the Hamiltonian matrix 
elements between physical states should be finite. From these conditions, we may derive the renormalization for 
the couplings and the boundary conditions for the meson wavefunction. The coupling constants are renormalized as 
follows 

4tt 4tt 1 f 1 1 4tt 4tt 1 f 1 1 

a?N ~ opV ~2J x(T^x)' a[N ~ ~ 2 J x(l - x) { ' 

As in the Gross-Neveu model, these integrals have been regularized which is not explicitly denoted here. It should be 
noted here that this generalizes the renormalization of the coupling constant in the Gross-Neveu model The mass 
parameter M needs no renormalization. This is again consistent with the renormalization in both the Gross-Neveu 
model and the 't Hooft model. We expand the meson wavefunction as 

<p(x) = <p {0) + tp {1) (l - 2x) + ip(x) (29) 

where (p^ 1 * are constants and <p(x)/ [x(l — x)] is integrable at x = 0, 1. Then, the boundary conditions for the 
meson wavefunction are 




b\ 1 
b+ \ 1 - 2x 



(30) 



Here, we defined b± = {\fP\ ± /2. In particular, when the coupling constants are equal, a 2 = a\ : or when the 
masses are equal, (3\ = P2, the boundary conditions simplify to 
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The meson wavefunction does not vanish at the boundaries. This property is similar to that of the Gross-Neveu model 
but unlike that of the 't Hooft model. When the Gross-Neveu couplings are zero, the wavefunction does vanish at the 
boundaries, thereby recovering the boundary conditions of 't Hooft. Also, it is instructive to check that the equation 
for the meson wavefunction (|2(]) and the boundary conditions (|ll) for /3i = (5i reduce exactly to the equations (0) in 
the Gross-Neveu model for a and \ scalars when tp^ = and t/^ 1 ) = 0, respectively. 
The equation for the meson states is now reduced to 

Hip(x) = n 2 f(x) 

- + ^) Viz) - P fdy + 2^ f- A + P 2 + In 1Z£) (32) 

x l-xj J [y-xr \ x J 

subject to the boundary conditions (^). It is straightforward to check that the Hamiltonian is Hermitean under the 
given boundary condition. The explicit dependence on the coupling constants is contained in the non-trivial boundary 
conditions. The problem has been reduced to that of solving a well defined integral equation. Below, we restrict to 
the case of scalar and pseudo scalar couplings being equal, a = 0,5, for simplicity. We will, however, consider fermions 
of different masses in general. The more general case may be dealt with in a similar fashion, involving somewhat more 
complicated formulas. From here on, we shall use the notation G = a 2 N/{A-K) to avoid cluttering the formulas. For 
the case a = 05, the matrix elements of the Hamiltonian may be simplified to the following form which is appropriate 
for the application of variational methods. 
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IV. SYSTEMATIC METHODS FOR SOLVING THE MESON BOUND STATE EQUATION 

In this section, we show how the meson bound state equation (|2^) may be solved systematically utilizing a finite 
dimensional system of algebraic equations. We will give explicit formulas for two approaches, namely a variational 
method using polynomials of the momentum fraction and a method of working with the equation directly using a 
sinusoidal basis. These methods will be used in the next section to investigate some physical properties of the gauged 
four fermi model. 



A. Variational method 

Here, we shall use a variational method using polynomials of the momentum fraction x that satisfy the boundary 
condition (p0|). In the 't Hooft model, a similar method was employed in ||. Without loss of generality, the basis 
may be chosen to be 



¥>2k(x) = G + 



B(k,k) 



¥>2k+i(x) = (1 - 2x) 



G (2k + l)[x(l-x)Y 
1 + 4G + B(k,k) 



fc = l,2,. 



(34) 



The normalization factor was chosen so as to make the matrix elements be of order one. The meson bound state 
equation becomes 



det{n 2 N kl -H k i) -0, H kl = (ip k ,Hifi), N k i = (cp k ,cpi) fc,/ = 2,3,4. 



(35) 



In practice, we approximate the solution by using a finite dimensional version of this equation. The matrix elements 
can be computed to be 



H2k,2l 
H2k+l,2l + l 



l fo a n sB(k + l,k + l) G a . kl 
2 + ft - 2) B(k, k)B(l, I) + 2 ^ + ^ + WTT) 
1 _ {2k + l){2l + l)B(k + l,k + l) G 

2 [Pl '' 2 ' (2k + 2l + l)B(k,k)B(l,l) 2(1 + 4G) 



(ft + ft) 



G 



1 + 4G 



G 



k 



I 



1 



H2k,2l + l — #2i+l,2& — ^ (ft _ ft) 



kl(2k+ l)(2l + 1) 
\ +iG\k + l ' l + lj 2(k + l)(k + l + l) 
G 1 (2l+l)B{k + l,k + l) 



l + 4G2fc + l (2k + 2l + l)B(k,k)B(l,l) 



(36) 



2fc,2Z 



(k + l)B(k + l,k + l) 



G 2 + G 



2(2fc + l) 2(2/ + 1) 



N- 



1 



2fc+l,2i+l — 



G 



2(2fc + 2/ + l)5(fc, &)£(/,/) 

(2fc + 1)(2Z + l)(fc + l)B(k + l,k + l) 
2(2k + 21 + 3)(2fc + 2/ + l)5(fc, £;)£(/, I) ' 3 Vl + 4G 
G ( k I 



2(1 + 4G) V2fc + 3 21 + 3 

N 2kl 2l+1 = N 2k +l,2l = 



(37) 



When the fermion masses are equal, (3\ = fa, the even and the odd sectors decouple. 
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B. Multhopp's method 



Rather than using a variational method, we may choose to work with the equation (|32|) directly. By a clever 
choice of basis, this singular integral equation may be brought into an algebraic equation. The method we use here 
generalizes the methods used to numerically analyze the 't Hooft equation previously jl0],f7| . 

We expand the meson wavefunction as 

<p(x) = ipW - tpW cos6 + ^2 Cn sin n6, x = i(l + cos 9) (38) 

^(0,1) terms are absent in the 't Hooft equation and can be determined in our model from the boundary conditions 
©• 

p<°> = 2nG ]T Cn, ^ = -27T— ^— c - (39) 



1 + 4G ^ 

n: odd n: even 



then the meson equation (B2J) becomes 



E 



[ti 2 P n (0)-M n (6)\ c„ =0 (40) 



where 



. G n: odd 

P n {6) EEsinne» + 27T<^ g 



i+4G cos " n: even 



- ,„x / A-l &-1 \ „ nsinn(9 f0 n: odd 

^ 2 (iT^ + 1^ J + ^"UST + i ^ (A " A + ^ S) - even W 

In practice, we truncate the basis to a finite dimensional one and systematically analyze the convergence of the 
solution by varying the dimension of this finite space, which we shall denote by K. Then the meson bound state 
equation may be reduced to a generalized eigenvalue problem. 

K 



^ [Ll 2 Pmn ~ M mn ] C n = (42) 
n=l 

where 

K K 

In this work, we choose g m {8) = 2 sin md/(K + 1) as in the 't Hooft model. Other choices may be more convenient 
depending on the parameters of the model. Similarly to the variational method given above, M mn , P mn oc (fli — (3 2 ) 
when m — n = 1 (mod 2) and the even and odd sectors decouple completely when (3\ = (3% . This property is preserved 
for any choice of g m (&i) as long as the property g m (&i) = (—l) m+1 g m (9K+i-i) is preserved. Unlike the variational 
method, however, the approximate solution obtained by truncating to finite dimensional solution space needs not be 
an upper bound on the true solution and in general will not be. 



V. PHYSICS OF THE GAUGED FOUR FERMI MODEL 



First, we would like to determine the parameter region where the behavior of the system is physically reasonable. 
In particular, the meson bound state should not be tachyonic. Here, we will perform the analysis for fermions of 
equal mass since we expect tachyonic mesons in the flavor singlet channel if tachyonic mesons exist at all. Using the 
variational method for the meson wavefunction in a manner similar to the previous section, 

Vl {x) = G + 5( 7 , 7)" 1 [x{l - x)] 1 7 >0 (44) 
we obtain an upper bound for the meson mass squared for each 7 > 0. 
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M 7 



, 5(27,27) 
B(7,7) 2 



(<Pi,<fh) G 2 + G^- T + 




(45) 



27+1 ~ 47+1 -8(7,7) 



This immediately establishes that when (3G < tachyonic mesons will appear. While it is logically possible from this 
analysis that the region with both /3 < and G < may be physically consistent, it is unreasonable to expect so; 
in practice, we find that tachyonic mesons appear for this case also, when we have a large enough variational space. 
Therefore, we henceforth interest ourselves in the region G > and j3 > 0. 

Using the methods explained in the previous section, we may obtain the spectrum and the meson wavefunctions. 
We plot the spectrum and the wave functions for some typical cases below in figures fig. |(| fig. and fig. |^. We have 
added a brief note on the convergence of the numerical data as an appendix. 
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FIG. 6. Mass squared of the lightest four meson states versus the coupling G for j3 = 1. 
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e- 




X X 
FIG. 8. The meson wavefunctions for the gauged four fermi model (solid), the 't Hooft model (long dashes) and the 
Gross-Neveu model (short dashes). The functions are plotted for the lightest two meson states. The parameters chosen 
for the gauged four fermi model is G = l,/3 = 1 for both mesons. For the 't Hooft model (3 = 1 also. For the Gross-Neveu 
model, in the lightest meson case, the meson mass was chosen to agree with that of the gauged four fermi model. In the next 
lightest meson case, fJ? / /3 = 3 was chosen. 



As in the 't Hooft model, the fermions are confined and the following Regge-type behavior is observed for the highly 
massive mesons, similarly to the 't Hooft model: 

- 7r 2 fc k > 1, (46) 

This Regge behavior is shown in fig. ^ 
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15 20 25 30 35 40 



FIG. 9. Mass squared normalized by 7r 2 of the meson states versus the "level number" for the case /3 = 1, G = 1. □ (+)'s 
denote states whose wavefunctions are even (odd) under x <-+ 1 — x, 



We may understand the behavior of the spectrum in the various limits of the model as follows: when we turn off 
the Gross-Neveu coupling a 2 , the spectrum reduces to that of the 't Hooft model. As we take the gauge coupling 
to zero, which effectively is the limit f3 — > oo, the splitting between the higher levels disappear. We have explicitly 
checked that the mass of the lightest meson approaches to the Gross-Neveu model value plotted in fig. [| in this limit. 
For the higher levels, \j? approaches 4/3 in this Gross-Neveu limit. The chiral limit may be identified as the limit 
G (or a 2 N) — > oo, similarly to the Gross-Neveu model case. In the spectrum, mass of the meson states decrease as 
we approach the chiral limit and it is clear that — s- for the lightest meson bound state in the limit G — > oo. When 
the limit (3 — > oo is taken in addition, it may be explicitly checked that the next lightest meson satisfies fj, 2 /M 2 — > 4 
corresponding the the a mass in the Gross-Neveu model. 

As (3 becomes large, the meson masses behave linearly with the quark masses, as is expected from the naive 
constituent quark picture. This picture is supposedly valid for highly massive quarks. The lightest meson behaves 
in a qualitatively different manner from the other meson states in the theory. This is a feature of the gauged four 
fermi model; such a behavior does not occur in the 't Hooft model. The Gross-Neveu coupling affects the lightest 
meson state relatively more than the other meson states. This disparate behavior is a necessary consequence of the 
Gross-Neveu limit where fi 2 / ' (3 of x and the other mesons approach the corresponding value fl 2 in the Gross-Neveu 
model and four. In the 't Hooft model, ^ 2 / (3 for the lightest meson also approaches four for large (3. 



VI. SUMMARY AND DISCUSSIONS 



In this work, we have solved the meson sector of the gauged four fermi model in the large N limit. We provided 
detailed prescriptions for solving the model systematically which should be useful for further work involving this class 
of models. We also determined the physically consistent region in the gauged four fermi model and analyzed the 
model there. However, it is possible that the regions with tachyons correspond to other phases of the model that is 
inaccessible to our current methods. 

Both the Gross-Neveu model and the 't Hooft model have been used extensively in the literature to understand 
the physical behavior of real systems, such as QCD. A model that combines the two models should be quite useful 
for understanding the dynamics of various field theories. In one direction, the four fermi coupling has been used to 
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model strong interaction dynamics involving dynamical symmetry breaking for quite some time |ll| ] . By interpolating 
between these two models, we intend to shed more light on the relation between the physical behavior of these two 
theories. Furthermore, when dynamical symmetry breaking scales are widely separated, in the intermediate energy 
range, the theory effectively becomes a gauged four fermi model. Such situations can occur quite generally where the 
lower energy scale is the QCD scale or some technicolor scale. These types of theories are of phcnomenological interest 
and have been studied actively, for instance, in the so called "top quark condensation" models jl2],[ll]]. Admittedly, 
the model we studied is a (1 + 1) dimensional toy model version of such theories. Historically, however, (1 + 1) 
dimensional theories have played an important role in understanding of the corresponding higher dimensional theories 
and we hope this will also be true in the future. 

Acknowledgments: We would like to thank Tomoyasu Ichihara for his collaboration during the early stages of this 
work. We would also like to thank K. Itakura and H. Sonoda for discussions. 



Appendix: A brief note on the convergence of numerical solutions 

The convergence of the numerical solutions depends on the parameters (G,f3). When j3 > 1, it is relatively easy to 
achieve a relative accuracy of ~ 10~ 4 in the meson mass at least using both the variational method (dimension ~ 10) 
and Multhopp's method (dimension ~ 400). For (3 < 1, more effort is needed to achieve the sam e leve l of convergence. 
The difficulties in the variational method using polynomials of the momentum fraction (section IV A) arise because of 
the round off errors since the eigenvalues in the normalization matrix (|37| ) tend to become small. Analytically choosing 
an orthonormal basis or using some other basis appropriate for the parameter region in question might alleviate this 
problem. In Multhopp's method (section |lV B| ), the limitations arise due to the necessary computational time when 
using larger space of functions. Choice of g m (Q) may speed up the convergence process in some parameter regions. 

In the parameter regions where the convergence is slow, extrapolation in the data can be effective. We have found 
that trial functions of the type x + aK b fits the data quite well. Here x is the extrapolated value, K is the dimension 
of the space span by the basis and a, b are parameters. Extrapolation can sometimes be misleading so checks on 
the results are desirable. In our case, we compare the extrapolation values from both the variational method and 
Multhopp's method and we confirm that they are consistent within the errors of the fit. An example of such an 
extrapolation is shown in fig. fTq. 
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FIG. 10. Extrapolation of the numerical data for the variational method and Multhopp's method. In this example, G = 1 
and (5 = 0.1 and = 0.095 ± 0.002. 
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